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A coupled model of magnetic flux generation and transport in stars* 
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We present a combined model for magnetic field generation and transport in cool stars with outer convection zones. The 
mean toroidal magnetic field, which is generated by a cyclic thin-layer aQ, dynamo at the bottom of the convection zone is 
taken to determine the emergence probability of magnetic flux tubes in the photosphere. Following the nonlinear rise of the 
unstable thin flux tubes, emergence latitudes and tilt angles of bipolar magnetic regions are determined. These quantities 
are put into a surface flux transport model, which simulates the surface evolution of magnetic flux under the effects of 
large-scale flows and turbulent diffusion. First results are discussed for the case of the Sun and for more rapidly rotating 
solar- type stars. 
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1 Introduction 

The increasing observational knowledge on stellar mag- 
netic activity makes it possible to use the observed ac- 
tivity patterns to constrain stellar dynamo models (Strass- 
meier 120051 Collier Cameron [20071 Holzwarth, Mackay & 
Jardine 2007). On the one hand, many solar and stellar dy- 
namo models in the literature make the assumption that the 
toroidal magnetic fields created in stellar interiors repre- 
sent also the surface emergence patterns. On the other hand, 
numerical simulations of the rise of flux tubes in the con- 
vection zone (Caligari, Moreno-Insertis, & Schussler 1995 ) 
have been successful in reproducing the observed prop- 
erties of sunspot groups, among which e.g., the tilt an- 
gle is of particular importance in surface flux transport 
models and Babcock-Leighton-type dynamos (Dikpati & 
Charbonneau 1 1 999b . Therefore, it is now possible to con- 
sider in a consistent way the interrelations between the dy- 
namo mechanism operating in stellar interiors, transport 
of toroidal magnetic flux in the convection zone, and the 
emerging flux which evolves under the effects of surface 
flows. Our approach combines models for three processes: 
1) the dynamo operating in the overshoot layer at the bottom 
of the convection zone, 2) stability and rise of magnetic flux 
tubes through the convection zone, 3) transport of magnetic 
flux on the surface. 



2 The model 

We consider stars with solar internal structure and differen- 
tial rotation, keeping the difference between the minimum 
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and maximum angular velocities, A £7, independent of the 
stellar rotation rate, which is a plausible assumption for 
rapidly rotating stars (e.g., Barnes et al. 12005b . The only 
stellar parameter that we vary is the surface rotation period 
at the equator, P rot . 

The dynamo model. The generation of magnetic flux 
is described by a dynamo model, which provides the 
toroidal component of the mean magnetic field. As a sim- 
ple example, here we consider an aft dynamo operat- 
ing in a thin layer at the bottom of the convection zone 
(Schmitt & Schussler 11989b . We assume a Sun-like radial 
differential rotation profile, according to helioseismic obser- 
vations (Schou et al. 11998b . A negative a-effect is assumed 
for latitudes below 35°, where the radial shear is positive. 
This choice is motivated by an a-effect due to flux tube in- 
stabilities (Ferriz Mas, Schmitt & Schussler 1 1994b or unsta- 
ble magnetos trophic waves (Schmitt 2003). The maximum 
strength of the a-effect (at A=±17.5°) is assumed to be pro- 
portional to the stellar equatorial rotation rate, Qq. The tur- 
bulent diffusivity of 5.6 • 10 11 cm 2 s _1 is chosen such that 
the cycle period becomes 1 1 years in the case of the Sun. 

The rise of flux tubes in the convection zone. In this 
part of the model, the thin flux tube approximation is con- 
sidered. Magnetic flux tubes are assumed to form (e.g., 
by magnetic Rayleigh-Taylor instability) within a layer of 
toroidal magnetic field created by the radial shear in the up- 
per tachocline. When the field strength of a tube exceeds a 
threshold, the Parker (undulatory) instability sets in, and the 
flux tube forms one or two rising loops. The probability for a 
flux tube to erupt from a given latitude A and at a given time 
t is assumed to be proportional to the mean toroidal field, 
B(\,t), provided by the dynamo model. The total number 
of erupting flux tubes per activity cycle is scaled with Qq, in 
accordance with the observed relation between rotation and 
activity (e.g., Montesinos et al. 12001b . Considering the lin- 
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Dynamo model vs. emergence pattern 




10 20 30 40 

time (yr) 

Fig. 1 Time-latitude diagram of the dynamo-generated 
mean toroidal magnetic field B(\,t) at the bottom of the 
convection zone (contours) and of the flux loops emerging at 
the surface (dots) for a Sun-like star with P rot = 27 d. The 
dynamo waves and the emergence pattern closely match in 
this case. 

ear stability of flux tubes in the mid-overshoot region (Ferriz 
Mas & Schussler 1 19951) , unstable flux tubes near the stabil- 
ity limit are chosen. The corresponding field strengths are 
around 10 5 G for the case of the Sun. Having determined 
the eruption times, initial latitudes, and the corresponding 
field strengths of the individual flux tubes, we perform sim- 
ulations of their buoyant rise (cf. Caligari et al. 1995} in or- 
der to determine the emergence latitudes and the tilt angles 
of the resulting bipolar magnetic regions at the surface. The 
separation of large bipolar regions on the Sun corresponds 
to much higher azimuthal wavenumbers (m = 10 — 60) than 
found for Parker-unstable flux tubes (m = 1 — 2). A solution 
to this problem is not available yet. There are, however, in- 
dications that bipolar regions dynamically disconnect from 
their roots at depths of only a few Mm below the surface and 
at an early stage of emergence, and this limits the azimuthal 
separation of emerging bipoles (Schussler & Rempel 120051) . 

Surface flux transport. The emergence times, latitudes, 
and tilt angles resulting from the previous step are taken as 
the input (source term) for a surface flux transport model 
(Baumann et al. 2004, 2006). The source term defines bipo- 
lar magnetic regions as a function of position and time. The 
flux transport code simulates the evolution of purely vertical 
magnetic fields at the surface, under the effects of Sun-like 
latitudinal differential rotation, meridional flow, and super- 
granular diffusion. The emergence longitudes of bipolar re- 
gions are assumed to be randomly distributed. The areas of 
the bipolar regions are determined at random with a number 
distribution N(A) ~ A~ 2 , which represents the observed 
distribution of solar bipolar magnetic regions (Schrijver & 
Harvev fT994b . 

3 Results 

First we consider the Sun-like case with P rot = 27 d. 
Figure Q] shows the time-latitude diagram of the dynamo- 
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Fig. 2 (a) Time-latitude diagram of the azimuthally aver- 
aged strength of the radial surface magnetic field for the so- 
lar model, (b) Time variation of the surface-integrated mag- 
netic flux. The values are averaged over 27-day time inter- 
vals. 



generated toroidal magnetic field at the bottom of the con- 
vection zone and the locations of flux emergence. The 
toroidal field contours almost coincide with the overall 
emergence pattern, because of the small poleward deflection 
of rising flux tubes for a slow rotator like the Sun. This result 
justifies an implicit assumption that is often made when in- 
terpreting solar dynamo models: the surface activity pattern 
reflects the dynamo wave pattern. Azimuthal averages of the 
radial surface magnetic field are shown in a time-latitude di- 
agram in Fig. [2k and the time evolution of the total surface 
magnetic flux is given in Fig. [3). The values are close to the 
results of Baumann et al. ([2004), with the exception that po- 
lar fields are slightly weaker in our case. The reason is that 
the tilt angles resulting from the flux tube rise are system- 
atically smaller than for the functional latitude dependence 
assumed by Baumann et al. (12004b . but actually provide a 
better match to the observations (see Fig. 12 in Caligari et 
al. 11995b . An animation of the evolving surface flux is given 
in the supplementary file Bsurf_2 7d.gif. 

We next consider a star with a rotation period of 10 d. 
The generated toroidal field and the emergence patterns are 
shown in Fig. [3 In this case, the poleward deflection of ris- 
ing flux tubes is stronger than in the case of the Sun. The 
responsible effect is the component of the Coriolis force, di- 
rected towards the rotation axis (resulting from a flow along 
the flux tube in order to conserve angular momentum), and 
is proportional to the rotation rate. As shown in Fig.[4k, we 
have strong polar fields (note the different saturation level 
of the colour table, as compared to Fig. [2^). This is due to 
1) higher emergence rate, 2) larger tilt angles owing to faster 
rotation. An animation of the evolving surface flux is given 
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Fig. 3 Same as Fig. [Q for P rot = 10 d. 
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Fig. 4 Same as Fig. for P rot = 10 d. 

in the file Bsurf _10d . gif . Figure |4j) shows that a cycle 
signal is no longer clearly observable when we consider the 
evolving surface magnetic flux. The reasons are 1) a larger 
degree of overlap between consecutive cycles owing to a 
stronger a-effect, and 2) strong polar magnetic fields in an- 
tiphase with the emerging flux. 

When we decrease the rotation period to 2 d, the pole- 
ward deflection of flux tubes causes a significant difference 
between the dynamo waves and the surface emergence pat- 
tern, as shown in Fig. \5\ An animation of surface trans- 
port is given in the supplementary file Bsurf_2d.gif. 
The tilt angles of emerging bipolar regions are much larger 
(around 35°) than in the solar case, leading to two latitudinal 
belts of opposite magnetic polarity. Meridional transport of 
the high-latitude belt leads to strong polar fields, with field 
strengths up to 30 G. 

4 Conclusions and outlook 

We have developed a consistently coupled model of mag- 
netic field generation and transport in the Sun and other 




Fig. 5 Same as Fig.(3] for P rot = 2 d. The dots are shown 
with one-month intervals. The emergence pattern is com- 
pletely different in comparison with the dynamo wave. 



cool stars. Similar to Schrijver & Title (2001) and I§ik, 
Schussler & Solanki (12007b , we find significant polar mag- 
netic regions for rotation periods of 10 d and 2 d, in addition 
to low-latitude activity. For P rot = 10 d, the surface flux 
transport blurs the periodic signal from the dynamo model. 
This indicates that for some stars the cyclic dynamo may 
be hidden in a non-cyclic surface activity, which is led by 
the combined effects of the dynamo (cycle overlap), large 
tilt angles, and the meridional transport. Furthermore, for 
P rot = 2 d we have found that the dynamo wave pattern 
and the emerging surface flux are completely different, ow- 
ing to the strong poleward deflection of rising flux tubes by 
the Coriolis force. 

As next steps, we plan to quantify the loss of magnetic 
flux from the dynamo layer in a consistent manner and con- 
sider it in the dynamo equations as a nonlinear term. Fur- 
thermore, we intend to develop a scheme which includes 
the evolving surface flux into a two dimensional flux trans- 
port dynamo as the source of the poloidal component at the 
bottom of the convection zone. 
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